EXISTENCE AND CONTINUOUS APPROXIMATION 
OF SMALL AMPLITUDE BREATHERS 
IN ID AND 2D KLEIN GORDON LATTICES 

DARIO BAMBUSI, SIMONE PALE AM, AND TIZIANO PENATI 

Abstract. We construct small amplitude breathers in ID and 2D Klein-Gordon infinite lat- 
tices. We also show that the breathers are well approximated by the ground state of the 
nonlinear Schrodinger equation. The result is obtained by exploiting the relation between the 
Klein Gordon lattice and the discrete Non Linear Schrodinger lattice. The proof is based on a 
Lyapunov-Schmidt decomposition and continuum approximation techniques introduced in [9], 
actually using its main result as an important lemma. 



1. Introduction 

In the present paper we prove existence of small amplitude breathers in some nonlinear Klein 
Gordon Lattices in dimension one and two. Moreover we prove that such solutions are well 
approximated by the ground state of a suitable nonlinear Schrodinger equation. 

The present paper is a direct continuation of [9] where the same result was obtained for the 
discrete nonlinear Schrodinger equation. More precisely, in [9] two of us proved the existence of spa- 
tially localized, time periodic solutions in one and two dimensional discrete Nonlinear Schrodinger 
equation (dNLS). In particular the so called Sievers-Takeno (ST) and Page (P) modes in ID, and 
also the hybrid (H) modes in 2D were constructed. The breathers of [9] have been obtained as 
critical points of the Energy functional constrained to the surface of constant I 2 norm. In turn they 
were constructed by continuation from the ground state of the Nonlinear Schrodinger equations 
(NLS). Thus such breathers turned out to be well approximated by the corresponding solutions 
of the continuous model. A key role in the proof was played by the interpolation of sequences 
(configurations of the discrete system) by the use of the so called Finite Elements, usually used in 
numerical analysis. 

In order to apply those ideas to the Klein Gordon lattice in which the trivial variational char- 
acterization of the breathers is no more true (see however [28] for general results on periodic 
orbits via variational methods, and [2,3] for recent results on KG chains), we have to establish a 
new connection between the KG lattice and the dNLS lattice. Actually the dNLS is heuristically 
known to be a resonant normal form (or modulation equation) of the KG lattice, an idea which 
has been exploited in order to describe the finite time dynamics of KG lattices. However, we 
need here a connection suitable for the description of infinite time dynamics. Such a connection is 
obtained by using the method of Lyapunov-Schmidt decomposition, in which the so called Kernel 
equation turns out to be a perturbation of the dNLS. We recall that a connection between normal 
form theory and the method of Lyapunov-Schmidt decomposition was first recognized in [13] and 
exploited in [7,8]. 

From a technical point of view the procedure is quite delicate, and in order to obtain a mean- 
ingful result we have to exploit the techniques introduced in [9] . We recall that in such a paper the 
idea was to use the Finite Elements in order to interpolate a sequence (configuration of the lattice) 
with a function of class H 1 . This allows to consider functionals on the discrete configuration space 
as restrictions to suitable subspaces of functionals on the continuous phase space. 
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This paper is part of a research line consisting in studying the dynamics of lattices using 
the continuous approximation. Along this line many results on the finite time dynamics have 
been proved [5,6,10,25,26,30], but little is known on the approximation for infinite times (see 
[14,15,17,29] on the dynamics close to solitons in FPU type models, and the papers [9,18-24,32] 
for what concerns existence of breathers). 

We recall that existence of breathers in lattices has been proved in [27] (see also [4,31]), by 
looking at the opposite limit, namely the anti-continuous one which leads to large amplitude 
breathers. Concerning small amplitude breathers, their existence have been established in one 
dimensional lattices using the spatial dynamics approach (see, e.g., [18-24,32]). In the case of 
higher dimensional lattices we recall the result by Aubry, Kopidakis, Kadelburg [1] (which in 
principle should be applicable also to the present model) and the results by Weinstein [33] (not 
directly applicable to the present model). However all these results are obtained by topological 
methods and give no information on the number of existing breathers and on their shape. The 
result of the present paper is actually the first one allowing to explicitly construct the ST, as well 
as the P and the H modes. 



2. Main result 

We consider the equation: 
(1) q J =a(Aq) 3 -V'(q 3 ) , j G Z™ 

with 

ns) = T"2^T2 |s|2p+2 ' {Aq)r - = £ 

P k:\k-j\ = l 

with |fc| = Em=l l fc m|< 

Here, and in the following, we will always consider n <E {1,2}. 

Remark 2.1. The value of the constant (3 is not essential since it can be changed by means of an 
amplitude rescaling. For our convenience we fix it as j3 = cj -1 , with ci defined in (|18[) . 



Remark 2.2. Due to the choice of the sign in front of the nonlinearity we obtain that the continuous 
approximation of the normal form is the focusing NLS equation. This is crucial for our analysis 
which does not apply to the defocusing case. 

To state the approximation part of our result, we need to refer to the ground state of the 
Nonlinear Schrodinger equation. Consider the Nonlinear Schrodinger Equation (NLS) 

lip = -A C V> - \ip\ 2p ip 



in K™, where 



" d 2 



is the usual Laplacian operator. It is well known that, if p < 2/n there exists a unique ground 
state ip c of the NLS fulfilling the additional requirements of being real valued, positive, radially 
symmetric and exponentially decaying [11, 12, 16]. Such a ground state is defined as the function 

which realizes the minimum of J Rn |VV>| 2 — ^xl'0| 2p+2 restricted to J Rn \ip\ 2 = 1. It can often 
be computed or described quite explicitly. 

For any fi > small enough consider the following 2™ distinct sequences ip l = {V'JO'O} - gzn 
defined by restricting the NLS ground state tp c onto Z™ as follows 



m M:=lM«). n = 1 . 



'^h,h := V'c (MUM*) , 

y£j 2 : =V>e (m'i,MJ2 + §), 

^juh '■= V»c 0-y'i + f ,Mi) , 
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These reference sequences correspond to the ST and P modes, plus the H modes in two dimensions. 

As a last step, for all the previous sequences, we renormalize the amplitude and add a temporal 
dependence in the following way: 

(3) = fip cos(wi)V/ , 

with u! = := y/l — mfi 2 , where m is a real constant (see sect. 13. ip . 

Remark 2.3. By construction the sequences are uniformly bounded in fj,, but with diverging 
£ 2 norm as \i — > 0; so we are calling "breathers" solutions which are localized on an increasing 
interval [—k,k] n with k ~ l/fi. The reference solution © share the same localization property, 
but with bounded £ 2 norm, due to the amplitude rescaling. 

We are now ready to state our result. 

Theorem 2.4. Assume n G {1, 2}, < a < \ and h <P < —> then there exists u* > 0, such that 
for any < /i < u* there exist 2" distinct real valued sequences q l = {<?}} - eZ „ € H 2 ([0,T];£ 2 ), 
which are time periodic solutions of fl} with period T = -^j^ . Such solutions fulfill 

In 

(4) h i -n HHW ) ^ c ^ r > r := p-2 +1; 

1 n 3 

(5) sup|«}(t)-*}(t)| <C 2A1 S , a :=___ + _. 

Remark 2.5. The periodic orbits we find are actually C 3 ' 1 in time, thus they are classical solutions. 
Indeed, once we get they are H 2 , and consequently C 1 ' by Sobolev embeddings, since the operator 
q — > \q\ 2p q maps C 1,1 into itself, one has q <E C 1,1 from equation ([T]). 

Remark 2.6. The bound p < 2/n comes from an analogous bound for the existence and approx- 
imability of the ground state of the dNLS (see [9]). On the other hand, for the dNLS, Wcinstcin 
proved in [33] the non existence of small amplitude ground states for p > ^. 

Remark 2.7. The reference solution, when measured in H 2 ([0,T];£ 2 ) has norm of order ft?" 1 , 
so the estimate (Q} shows that in such a norm the distance between the actual solution and the 
reference solution is small compared to the size of the solution. On the contrary such an estimate 
gives no information on the distance between the single particle in the approximate and the true 
solution (sup norm). A better and relevant control is given by the estimate ([5]), which is obtained 
through the use of a discrete analogue of Sobolev embedding theorems (see Sect. [6]). 

Remark 2.8. We also stress that our approximation estimates, not only control the spatial profile, 
but also the time dependence. In particular, one easily gets from (|27|) and (|3T|) that the first 
harmonic gives the principal contribution, the others being small corrections of order fi a , with 

° ■= p - 2 + 2 - 

The rest of the paper is devoted to the proof of theorem 12. 4[ and is organized as follows. 
Section [3] contains the setting of the problem, with the description of the proof of the main 
Theorem. The range equation is dealt with in Section [U while the solution of the kernel one is 
discussed in Section final estimates are presented in Section [5J Some technical details are given 
in Appendixes: in Appendix IA1 we show some regularity results for the nonlinearity; in Appendix IB1 
we give some improved estimates on the approximation of £ q norms; and in Appendix [C] we prove 
the extension of the Implicit Function Theorem used to prove Proposition [53] 

3. Settings and proof of Theorem 12.41 

In all the paper we will deal only with sequences qj which are reflection invariant, and thus 
which fulfill qj = q_j. Thus, when writing £ 2 we will actually mean the subspace of £ 2 composed 
by symmetric sequences. The same will be true for all the other spaces of sequences that we will 
meet in the paper. 
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More precisely, we will denote the space by £ 2 whenever it is endowed with its standard scalar 
product and norm, namely 



and it will be denoted by Q when endowed with the norm 

(6) Ikllg := \\qfi2 + -^(q, -Aq)j 2 , 

which will play a fundamental role in Lemma l4.5l and Section [SJ here \i > is the small parameter 
which was introduced in ([3]). 

Remark 3.1. From the technical viewpoint, the interplay between the norms ||-||^2 and ||-||g is 
one of the delicate points. Indeed it turns out that the Q norm is too strong to ensure enough 
regularity for all the continuation procedures; but the I 2 one is too week to grasp all the relevant 
information, in particular the non-degeneracy of the solutions, contained in the result of [9]. We 
will thus play with both the norms depending on the situations. 

We look for T-periodic solutions of ([1]) of the form 

qj(t) = Uj(u;t) , 

witfQ 

(7) u j(t) = ^ u j,i cos(it) . 

Wc define Xu C H k ([Q, 27r];i' 2 ), to be the space of even functions of time taking value in I 2 
which are square integrable together with their weak derivatives up to order k. 
In this paper we will only use k € {0, 2}; for k = 2 we will use the norm 

(8) \\uf X2 := (\Ht)\\% + \\u{t)\\l + \\m\i) dt = ]T(1 + I 2 + l A )u 2 , . 

Using the time-rescaled variable u, we rewrite ([T]) in the form 

(9) L^u = N(u) , 
where 

(10) : =L w -aA, := u? d tt + 1 ; 



are respectively the linear and nonlinear operators. 



(11) (N(u)) J : =(3\ Uj 



3.1. Lyapunov Schmidt decomposition. We look for small amplitude solutions, bifurcating 
from the trivial one. We perform a Lyapunov-Schmidt decomposition with respect to L\ = 
Lu\ v Define 

(12) V 2 := Ker{Li) , W := Range(L x ) , W 2 := W n X 2 

and Vo, the completion of V 2 in Xq. Such decomposition is invariant under the action of the linear 
operators (fTOj) . 

We correspondingly decompose u in the two components v and w 

(13) u = v + w, uEX 2 , veV 2 , weW 2 ; 

according to the development (J7J), if we denote e; = cos(lt) and C := {I S Z, I > 0, 1 ^ 1}, we may 
write 

(14) Uj(t) = VjCj + Wjjei , 

tec 



1 Due to the autonomous and reversible nature of ([TJ, it is rather natural to look for solutions even in time, 
thus with a Fourier development in cosine only. 
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where Vj = Vj x is the only Fourier component in the kernel. Denote A := /i 2 m = 1 — to 2 , ILy 
the projector onto Vq, and U\y — Ily — I the projector onto Wo- We remark that V2 and Vq are 
isometrically (up to a constant factor) isomorphic to I 2 , since any element v G V2 or z G Vq have 
only one Fourier component. 

3.2. Proof of Theorem 12.41 In this subsection we give the steps of the proof, leaving the full 
details to the subsequent sections. 

Step 0: Decomposition. We first decompose © in the following two equations, the first on the 
range and the second on the kernel 

(15) L^w = IL w (N(v + w)), 

(16) Xv-a(Av) = U v (N(v + w)). 

Step 1: Range equation. As usual in the Lyapunov-Schmidt decomposition, we first fix v and ui 
and solve the range equation (p~5|) via the Implicit Function Theorem (see Proposition ^. 2[) showing 
that 

w(v) ~ (l^J Hw(N(v)) + higher order terms. 

We stress that the range equation is solved using always the £ 2 norm. All the corresponding 
discussions and proofs are presented in Section 01 

Step 2: Kernel equation. We insert now the solution of the range equation w(v) in the kernel 
equation (JTSJ). Due to the smallness of w(v), the term HyN(v + w(v)) can be split into a main 
part TLyN(v) and a remainder: 



(17) Il v N(v + w(v)) = U v N(v) + Ily [N(v + w(v)) - N(v)) 
it turns out that (UvN(v(t)))j — (3ci\vj\ 2p Vj cos(i), and 

(18) cx := / \cos(t)\ 2p cos 2 {t)dt . 



By applying the following scaling 

(19) A = mfi 2 , v = (j,p<fi, 

factoring out the time dependence and recalling that we set j3 = the kernel equation then looks 

(20) -4(A0) + ™0-|^| 2 V + i?y(0) = O, 
where 

Since Ry (4>) is small, (see Lemma I5.3[) , the kernel equation, in the form (]20[) , appears as a 
perturbation of the dNLS model studied in paper [9] . The main result of [9] ensures the existence 
of breathers in the form of non-degenerate ground states. We exploit non-degeneracy to continue 
such solutions to solutions of the full equation (|20|) . 

In order to exploit the result of [9] we need to work using the norm Q rescaled by a factor 
so let us introduce the following notations 



(21) ^:=(>;l|-|Uj 



and 



(22) 2„:=(>;|H| Q J 



U :— (M n \\'\\fi , 



ls„^"IMI 2 s 



In terms of these norms it is possible to prove good estimates for Ry , see Lemma 15.31 In par- 
ticular we use them in order to prove some discrete analogue of the Sobolev embedding theorems. 
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This ensure the applicability of the Implicit Function Theorem I C . 1 1 for the continuation proce- 
dure (see Proposition 15 . 51 for the details); we thus obtain the 2™ solutions v l — fJ,p(f> 1 , with <\> % close 
to ip* (see estimate (|55|) and (14TH0 . and the corresponding 2™ solutions of (JTJ) 

q l {t) = v l (ujt) + wiv^ujt)) , i = 1, . . . , 2" . 

Step 3: Estimates. Finally one has to collect all the estimates in order to get the result. We will 
also obtain the improved estimate needed to control the sup norm of the difference between the 
true solution and the actual solution. This will be done in Sect. [5] since it actually requires the 
results of all the previous sections. □ 

4. The range equation 



In this Section we will prove Proposition 14.21 on the solutions of the range equation. We start 
by controlling the inverse of the linear operator L' u ' defined in (TH 



Lemma 4.1. If < a < h and \u 2 — 1| < \, the linear operator [lJ^\ is bounded from 
W C X to W 2 C X 2 

Proof. We use Neumann formula for the inversion of a linear operator. We rewrite as 

=L u (l + aL~ 1 A), 

from which 



(LM) =(I + aLZ 1 A)- 1 LZ 1 = J2(- aL ^ A ) k K X - 

k>0 

We observe that the series defines a bounded linear operator if 

( 23 ) a \K lA \\c(w ,w 2) < 1 - 

On one hand A acts only on the spatial index j and defines a bounded operator on £ 2 with 

|| Aw|| £2 < 2 \\w\\ e2 ; 

on the other hand, L^ 1 acts only on the temporal index I and, provided \u> 2 — 1| < 1/2, is bounded. 
Hence (|23|) is fulfilled provided a is small enough, i.e. a < h. 

□ 

We will show that w is smaller than v so that N(v + w) = A r (w)+higher order terms. It 
is thus natural to expect the solution of the range equation (fX5|) to be close to the solution of 
L^w = U w (N(v)) namely to 

w (v) := (L^ 1 U w N{v) . 

Define 

(24) U s ~{veV 2 : \\N(v)\y + \\N'(v)\\ c(e2iXo) <s} 

where we used the identification of V 2 with £ 2 (see (IT4)) ). 
We are now ready to prove the following proposition 

Proposition 4.2. Under the assumptions of Lemma \4- 1\ there exists 8 > 7 constants C, C , and 
a function w € C 1,:L (%, X 2 ) close to wq such that w — w(v) solves (|15l) . Moreover the following 
estimates hold 

(25) \\w(v)-Mv)\\x 2 < C\\N(v)\\ e2 \\N'(v)\\ 

(26) W(v)\\ c{X2>X2) < C'WN'^W^^y 

Remark 4.3. Clearly w(v) is small with N(v), i.e. 

(27) ||^)|| X2 <C"||Ar(t,)||, 2 . 
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Proof. Since the Nemitski operator defined by N is C 1 ' 1 (X 2 ,X ) by lemma lA.41 the implicit 
function theorem ensures the existence of a neighbourhood of the origin in which the function w is 
well defined. In order to get the size of such a neighbourhood and to prove the estimates (|25|) . ([26]) 
we go through the proof of the implicit function theorem using the contraction mapping principle. 
Let us first check that wq(v) ^ 0. We write explicitly wq(v): 



\n w N(v)]i = liw (M^jN^ex) = N 2 S"^- E [ r Iei| 2p eiezdt 

lec LJ ° 



= \ v j\ 2Pv j^Z^2 c l e l> 
lec 

since q is the Fourier coefficients of | cos(t)| 2p cos(f), there surely exists at least one q 7^ with 
I G C. 

Now rewrite the range equation as follows: 

(28) w=(l^A J1 w N(v + w) =: F(v + w) . 

Fix a positive Si, and v with \\v\\ < Si, and look for conditions such that, the r.h.s. of ([28]) defines 
a contraction of a ball of center wq(v) and radius S2 ■ We claim that 

\\w ~ w (v)\\ X2 < S 2 => \\F(v + w)-wo(v)\\ X2 <S 2 

if 82 ■— K ||iV(u)||^2 ||^ / (' u )||£(f2 x a ) an< ^ K sufficiently large. 
To prove it, first remark that, since N £ C 1 ' , 

(29) sup ||iV»|| £(X2 , Xo) < IIJV'MH^^,,) + C'S 2 , 

\\u-V-Wo(v)\\x 2 <$2 



from which 

\\F(v + w) -w (v)\\ x 



(l^A 1 n w [N(v + w) - N(v)} 
V ' x 2 

< C sup ||^V'(u)||c(a 2 ,Ao)I|w||a 2 < 

|m--u-io (-u)||x 2 < 5 2 

< C (C'S 2 + \\N'(v)\\ ci x 2 , Xo) ) (S 2 + C" \\N(v)\\ i2 ) = 

= ^ [C'K \\N{v)\\ l2 + 1) (k\\N'(v)\\ c{X2Xo) + C") S 2 , 

which is smaller than S 2 provided K is sufficiently big and ||iV(ti)|| £2 and ||AT'(w)||£^ 2 Xo ) suffi- 
ciently small, i.e. if v £ Us- Then one immediately sees that by possibly reducing S the r.h.s. of 
(f28l) has a Lipschitz constant smaller than one in the considered ball. So, we get the first of the 
estimates (|23|) . 

In order to get the estimate of the derivative of w just remark that 

w'(v) = (1 - F'(v + w(w)))" 1 F'(v + w(v)) 

where we denoted F(u) := [L^^j 1 HwN(u). Using Neumann formula in order to compute 
(1 — F'(v + ui(u))) -1 one sees that this is a well defined bounded linear operator provided 6 is 
small enough. Adding the estimate of F 1 in the ball, which in turn is obtained through (|2U)) one 
gets the thesis. □ 

Remark AA. Since the topology induced by the Q norm is stronger than the £ 2 topology, one also 
has that the solution w(v) of the range equation is C : (Q n Ug, X 2 ). 



In order to use in an effective way the inequality ([27]) we will make use of the following lemma, 
which will be proved in Appendix B, and which is based on the use of Sobolev embedding theorem 
applied to functionals interpolating the £ q norms. 
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Lemma 4.5. One has 

(so) m b EN 9 ^ c H« 



9 



which gives in the case of operator N 

(si) \\N(v)\\ el <c mii; 1 . 

Remark 4.6. The estimate (|31[) is much stronger than the trivial one obtained by using the ho- 
mogeneity of N and rescaling v. This will be crucial for our development. 

5. The kernel equation 
We have seen in formula (|20|) that the kernel equation looks like 
(32) G(<f>,ii) = 0, G(<M = G ! o(<M+iM<M 

where 



(33) G (<t>,fi) :=- — (A0)+m0-|0| 2 P0-O 



a 

is the equation for the ground state of the dNLS model. The maps G and Go will be considered 
as maps G : x R — > ^. 

The idea is to continue a solution of ([33]) to a solution of (|32|) . So, first we show that i?y is 
actually a perturbation of Go- Denoting by $ = $(/j.) a solution of Go = 0, we then show that 
G ($) is an isomorphism of Q M on 

We begin by recalling the statement of the main result of [9] . First we need to introduce a few 
objects. Consider the functional 



(34) Haffl := //' 

and the surface 



\j-k\=l\ j.kei n ^ jez™ 



5 = {0 : JV(0) = 1}, N{<j>) £ l& 

Theorem 5.1. (Main theorem of [9].) For any /i small enough and h < p < 2/n t/iere exist 2™ 
distinct real valued sequences $!(/x) which are solutions of (|20|) tuii/i i?y = 0. Smc/i solutions are 
coercive minima of Hq | „ . TTie coercivity is intended in the norm . Furthermore one has 

(35) ||$* — -0*11 2 ^ < C\x , 
where tpj are the sequences defined in ([2|). 

Remark 5.2. In the following we will concentrate on one of these solutions, so we will suppress 
the index i, from from \1/ and also from tp. 

Lemma 5.3. R v e C 1 ' 1 ^,^) fulfils 

(36) \\Rvm\q < UWX 1 , \\Rv{4>)\\c{q^d < c^uwX ■ 

Moreover, i?y(0,/i) and R' v ((j),(i) are continuous with respect to [i £ (0, +oo). 

Proof. The smoothness of Ry follows from the smoothness of N, since both the norms in the 
spaces Q and I 2 have been rescaled by the same factor. 

We prove the first of (|36| working first on the non rescaled quantity 



< 



\\U v [N(v + w(v))- N(v)}\\% < \\N{v + w{v))-N(v)\\ x < 



< C [\\VN(v)\\% + \\VN(w)\\%)\\w(v)f X2 



< 



C l (\\VN{v)\\% + \\VN(w)\\%) \\N(v)\\% 
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where we set ViV(u) := \v\ 2p . We now apply the scaling v = ^4>] and using Lemma [43] one has 



\\vn(v)\\% = ^Ki 4p <A n|1 ' r 



4p+2 
j 



= £N 4p+2 <^ 



If we ignore the term || ViV(io)||^ 2 which is much smaller in than the main one || ViV(u)||^ 2 and 
we take the scaled norm ||-L 2 we obtain 

\\Rv(4>)\\ei ^ c ^U\\t?- 

In order to conclude, we move to the estimate of R' v (v). Notice that 

\\ R vWh(Q^ei) = sup jjfcjf^ =sup jpijk H^vWIU((W a )- 

Let us first differentiate the non rescaled version of Ry- 

(37) — U v o [JV(u + w(v)) - iV(u)l = n v o [N'(v + w(v)) - N'(v)} + U v o N'(v + w(v)) o w'(v) , 
du 

We deal with the first addendum in the r.h.s. of (f3"T)l . by using the Lipschitz continuity of N' 

\\U V o [N'(v + «,(«)) - iV»]|| £(s/2) < ||[iV> + «;(«)) - iV'( V )]|| £(X2iXo) < 

< C'(IHIx a + IHIx a ) ap " 1 IHIx a < 

< CrdlHI^ + HiVWII^) 2 ^ 1 HiV^H^ . 
Hence, after rescaling the variable and the norm, we get 

c»\m P + f\\Nmi*? p - l \w{<t>)\\ii <c 1M v- np+ * <ci^w£ , 

where we have ignored the smaller term /i 2 ||iV(</>)||^2. A similar estimate can be obtained for the 
second addendum in (|37j) . Coming back to Ry 

the above informations yield 

The continuity with respect to the parameter fi follows from standard arguments. □ 

Lemma 5.4. Let $ fee a coercive minimum of Hq\ s , then the differential Gq($) is an isomorphism 
ofQ^oni*. 

Proof. By the theory of Lagrange multipliers one has that $ is a (free) critical point of Hq + mN, 
with a suitable m, while Gq($) is such that 

d 2 (H + mN)($)[h,h] = (& ($)h;h)e . 

Introduce now "polar coordinates" 

and write Gq(<I>) as a block matrix in terms of such coordinates. It has the structure 

'A b 
b T d ■ 

By non-degeneracy one has (A^Oi* > C||f||g , which, by Lax-Mil gram lemma implies that 
A : <2 M — > £ 2 is an isomorphism and is positive definite. 



(38) G' ($) 



in 
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Now one has 

(39) d = d 2 {H + mN) ($)[$,*] = (G'(#) [$],$) = -2p ^ ($) 2 p+ 2 <0 



which shows that the quadratic form d 2 (H + mJV) ($)[$, $] has a negative direction. It follows 
that such a quadratic form does not have null directions. Thus non-degeneracy and the thesis 
follow. □ 
Lemmas 15.31 and 15.41 give the hypothesis needed by Theorem IC.U We get the final 

Proposition 5.5. Let be a non-degenerate critical point of Ho\ s , then, for fj, small enough 

there exists a solution (f>(fi) of the rescaled kernel equation (|20[) . such that 

(40) u-n Q ^<c^. 

Proof. We verify the assumptions of Theorem lC.il Set X := and Z := £ 2 and xo(/x) := 
Define 

F :=G ((M, Fi:=R v (<t>,fi) . 

Then, by the previous lemmas the assumption of Theorem lC.il are satisfied and the thesis follows. 

□ 



6. End of the proof 
We begin a section by a simple lemma needed to obtain the estimate [5] 
Lemma 6.1. For any j G Z" we have 

(41) \4>i\<^H\\ a - 



Proof. We write the proof for the case n — 2. The case n = 1 is simpler. Denote j = (h, k); one 
has 

h h 

t>m,k ~ <Am-l,fc) = X (^ m < fe ~ t l>m-l,k)((t>m,k + <t>m-l,k) 
m— — oo m— — oo 



E 



which gives 



sup <p hk 
(h,fe)ez 2 



< 



< 



I 2 < 



m 6: 2 



r 2 



1/4N 



and E 



□ 

£W o/ £/ie proo/ o/ theorem \2.4\ We collect now the estimates needed to conclude the proof. From 
proposition 15.51 and theorem 15. II we get 

U-^PWq^ <c m . 

From this we get (using also Lemma 16.11 and Theorem 15. ip 



sup 

j 



v - ^ 1/J V 



\x 2 



i+2 



\wj\ < Cfi" 2 



from which the thesis immediately follows. 



□ 
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Appendix A. Smoothness of Nemitski operators. 

Lemma A.l. If p > then the operator N defined in (jTTJ) is C 1 ' 1 from I 2 to £ 2 with the usual 
norm. 

Proof. We first remark that I 2 C £ q for all q > 2, indeed 

HI!. - E N 9 = E(N 2 ) f ^ C ^J2 u i) f = c i IMI* > 

3 3 3 

which immediately tells that N(u) G I 2 if u G £ 2 . Indeed 

iiahii £ 2 2 = m 4p+2 = ni&s < c P iiuiir 2 . 

Moreover it will be useful to remind that £ 2 C t°° since 

IM| 2 oc. = (sup|Mj|) 2 = max|u,-| 2 < \u \ 2 = ||u|| 2 2 . 

The first continuous embedding immediately gives also the continuity at the origin. To obtain the 
continuity at a point « ^ we proceed showing that N is Frechet differentiable at any u with 
bounded differential N'(u). From a direct computations one has that 

N'{u)[h] = \u\ 2p h, 

hence 

\\N'(u)[h}\\% = M 4 ^ 2 < Ci |M| 4p \\h\\% 

3 

which gives 



\\N\u)\\ c[PP) <C p \\u\\ 2 v. 

The possibility of locally bounding the differential yields to the local Lipschitz continuity of N, 
since 

\\N(u) - N(v)\\ e2 < I sup l|A'HII £( ,2 /2) ] \\u- v\\ £2 . 

\\M\ e 2<\\u\\ i2 +\H\ e 2 7 

Finally, the (local) Lipschitz continuity of A' is due to p > \. Indeed 

\\N'(u)[h] N'(v)[h}\\% = J2 (N 2P - I^| 2P ) 2 h] <C p \\u- v\\% \\h\\% 

3 

with C p = C(p, Hull^oo , U^H^oo) a local constant. □ 

Corollary A.2. If p > \ then the operator N is C M (Q, £ 2 ). 

We are now interested in regularity of A as a map from X 2 to X . We first state an auxiliary 
Lemma: 

Lemma A.3. Let us define f(t) := \\u(t)\\ e2 : [0,2w] -> W . If u G X 2 then f G ^([O, 2tt]). More 
precisely one has 

(42) \\f\\ Hl <V2\\u\\ X2 . 

Proof. Surely / G L 2 (I), where we set X :— [0, 2tt]. By differentiating f 2 (t) we get 



thus 



\f(t)\ 2 <^(\\u\\ e2 \\u\\ e2 f <2\\u\\%; 



u 



£2 

hence 

/■27T /.27T 



\f{t)\ 2 dt<2 \\u\\l < 2||-' 1 







A' 2 
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which implies / £ The estimate (H2J) follows from 



11/11 



ff 1 



2tt 



2 71 



2tt 



\f(t)\ 2 dt + \f(t)\ 2 dt<2 \\u(t)\\t> + \\W)\\t> <2||«| 



IX, 



Lemma A. 4. TTie nonlinear operator N is C ' (X2, Xq). 

Proof. From Lemma [A. II it follows that, for any fixed t £ I it holds 

IliVKt))!!^^!!^)!!^ 2 , 

so from the Sobolev embedding L q (I) if 1 (J), g > 2 we get 



2,T 



2tt 



l|iv( u (t))||^ < Cx / |K*)lljr 2 < c 2 HHi^H^ < c 2 \\u\\ 



4p+2 



X 2 ' 



which simply gives 



\\N(u)\\ Xo <C\\u\\ 2 l +1 



Let us consider now the Frechet differential 

N'(u) :/ieI 2 H N'(u)[h] = \u\ 2p h e X ; 
once more, for any fixed t £ I one has 

\\N\umHt)]\\l<c 1 \\um\t\\h{t)\\% 

thus, again from Sobolev embeddings, we get 



2 77 



2tt 



< Ci / IK*))I# IIM*)I& < C2 \\u\\% \\h\\ X2 ; 
we have so proved that the differential is locally bounded 

II^MII £( x 2 ,x„)<C||«||^, 
and hence N £ C 1 (0). Gathering the previous results we deduce 

\\N(u)-N(v)\\ Xo <( sup \\N'(w)\\ c{X2!Xo) )\\u-v\\ X2 

\wG[u,v] J 

which is the local Lipschitz continuity. The local Lipschitz continuity of N'(u) can be obtained in 
the same way as in Lemma I A. II Indeed 

r2n 



\\(N'(v + w)-N'(v))[h}\\ 



< 
< 



\\h(\v+ w \ 2p -\v\ 2 n\\l< 

\\h\\%\\\v + w\ 2p -\v\ 2p \t < 



\\\v + w\^-\v\ 2p \ 



< 



L- 



< C\\h\\ 



X-2 



\v + w\ 2p - M 2p \ 



If 2 



L 2 



Moreover, following Lemma \A. II one has 



\\v + w\ 2p -\v\ 2p \\% <C L (tf\\w\\% 



with 



C L {t) = 2p{\\v{t)\\ toa + \\w{t)\\^fP- 1 < 

< 4(||||w|| <a || £00 - 

< 4(|H| Xa + HI^) 2p - 1 



w\\ l 4 L ~) 2p - 1 < 



So 



\\N'(v + w)-N'(v)\\ ciX2tXo) <L\\w\\ 



x 2 ) 



L:=C(\\v\\ X2 + \\w\\ x ) 



2p-l 
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where the Lipschitz constant L is local, unless for p = \. □ 



Appendix B. Approximation of discrete functionals and proof of Lemma |4~1j1 



The proof of Lemma [4.5l is based on the use of Sobolev embedding theorem applied to the contin- 
uous interpolation of some discrete functional. In turn, following [9] , the continuous interpolation 
is obtained through the method of finite elements as we are now going to recall. 

B.l. The case n = 1. Define the sequence of functions Sj(x) by 

0, if \x-j\>l 

(43) Sj(x) = { x-j + 1, if -l<x-j<0 

-x+j + 1, if < x ~ j < 1 

and, to a sequence ipj, we associate a function 

(44) T(x):=Y,4>j»j{xlti- 

i 

On the interval Tj :— [fij, Li{j + 1)) the above function reads 

(45) T( g ) = (x- M i) ( ^ +1 ~^ ) +V> J -. 

B.2. The case n = 2. For each multi-index j = (h,k), take the function Sh,k{x,y) which rep- 
resents the hexagonal pyramid of height one centered in (h, k) with support the union of the six 
triangles T^ k , T^_ lk , T^ k+1 , T^ k _ v T^ +lk . For example, on T^ k the function s h ^ k represents 
the plane in R 3 which passes through the three points (h, k, 1), (h + 1, k, 0), (h, k + 1, 0), namely 

s h<k (x, y) = -x - y + h + k + 1. 

We take {sh tk (x / 'fi,y / 'fx)} /h,k)e% 2 as a basis to generate a piecewise linear function ip(x,y) which 
interpolates iph,k 

(46) T(x,y) := ^ ip h ,kSh,k( x / ^ S//m)- 

(h,k)el? 

On the triangle Tj^ k the function ip is a plane which reads 

(47) T{x,y) = {x- fxh) h (y - fik) \-iph,k, 

H n 

while on the opposite triangle ' k it reads 

( AQ\ T ( \ I u A^h,k-^h-\,k) , M (lph,k ~ ll>h,k-l) . , 

(48) T(x,y) =(x- fih) V [y - fik) h iph,k- 

fi fi 

Definition B.l. We denote by £ M the linear subspace of H 1 of the functions (|46)l with {ipj} € Q^- 

B.3. Interpolation. We recall now some lemmas which were proved in [9]. 

Lemma B.2. Let Te^, denote by ip = {tpj} the corresponding sequence, then one has 

(49) / \VT{x)\d n x = -\fi n {iP,&^) P 

Jm. n 

Lemma B.3. Let T 6 E^, denote by tp = {ipj} the corresponding sequence; define 

G C (T) := [ |T|«+ 2 , G d (T) ■= a n V + 2 , R G (T) := G C (T) - G d (T); 

if q > 1 i/ien i?c G C 2 (£^) and for any bounded open set IA <Z there exists C{U) such that 

\\Rg\\ CH u) < 
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Proof of lemma \4-5\ One has 

(i n y2W 9+2 = I \T\ q+2 d n x + 0(fi) < C\\T\\*}? + 0{n) 

= c\\nZ + o^)<c\\nZ 

where the third inequality follows from the fact that, for n = 1,2 one has that the injection 
H 1 <—> L p is continuous for any p. From this the thesis immediately follows. □ 

Appendix C. A version of the implicit function theorem. 

Theorem C.l. Let X,Z be Banach spaces and let S I C EL Let F S C X,1 {U xI,Z), with 
U C X open. Let xq = Xo(fi) € U be such that 

(1) 

(50) F Q (x Q (ji),n) =0;V^eX 

(2) Fq x (x,/i) is Lipschitz inhi uniformly in fi, i.e. there exists L — L(U) independent of /i 
such that 

(51) \\F^ x (x,fi)-F^ x {xo,fi)\\ c{xz) <L\\x-xo\\ x , VxEU; 

(3) Fq x (xq, fi) is invertible and its inverse is bounded uniformly in fi, i.e. there exists C\{U) 
such that 

(52) WiF^ixo,^}- 1 ]^^ <Ci, V^el. 

(4) Let F\ 6 C l {U x X, Z) be such that there exist a > and Cz(U) such that 

(53) \\F 1 (x, t i)\\ z <C 2 n a , V(x,[i)eUxl; 

(5) there exist [3 > and C^{U) such that 

(54) \\Fl x (x,v)\\ c{xz) <C 3 ^, VMeMxI; 

Define F := Fq + F\ , then there exist Us C U and [J,* (6) < /i and a function x(n) : Xq := (0, /x*) i— * 
Us which solves 

(55) F(z(/i),/i)=0, 
wii/i 

(56) ||a:(/i) - ioMHjc 
Moreover one has x(fi) G C°(Xq,Us). 

Proof. First remark that, provided /i is small enough and possibly restricting U, [F^x)}^ 1 exists 
and fulfills 

(57) IW^rl^x) <C, V/xeJ,. 

Define now A(:r, /u) := x — [F x (xq, [j,)]~ 1 F(x, h) and remark that any fixed point of A is a solution 
of our problem. We now prove that A is a contraction of a ball of radius 0(fi a ) and center x$. 
So, let x be such that ||x — xq\\ < 5, and let us estimate the Lipschitz constant of A in such a ball. 
One has 

\\F^ x (x ,^-F^(x,ti)\\ <CS 
\\Fi tX (x o ,n)-F[Jx,fx)\\<Cn 

It follows that 

\\A'(x,ri\\ = \\[F x (x ^)}- 1 (F x (x ^)~F x (x^))\\<q<l 
provided /i and 5 are small enough. Compute now 

\\A(x,fi) - x \\ < \\A(x ,n) ~ x o\\ + \\M X ^) ~ M x o,n)\\ 

< Cfl a q6 , 
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which is smaller than 6 provided C/J, a < (1 — q)S, which in turn can be obtained e.g. by taking 

2(1 -g) ■ 

□ 
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